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Abstract 

We explicitly show the fully nonlinear equivalence of the <5N and the covariant 
formalisms for the superhorizon curvature perturbations, which enables us to safely 
evaluate the non-Gaussian quantities of the curvature perturbation in either for- 
malism. We also discuss isocurvature perturbations in the covariant formalism and 
clarify the relation between the fully nonlinear evolution of the curvature covector 
and that of the curvature perturbation for multiple interacting fluids. 



1 Introduction 



Inflation is now widely accepted as the mechanism of generating primordial density fluc- 
tuations, which are almost scale invariant, adiabatic, and Gaussian. Such fluctuations are 
generated inside the horizon and expanded to superhorizon scales due to the inflationary 
expansion of the Universe. These features are confirmed by the recent observations of 
cosmic microwave background anisotropies like WMAP experiments [I]. In particular, 
the large-scale anticorrelation of the temperature and the E-mode polarization strongly 
supports the presence of the superhorizon epochs of primordial curvature perturbations 
[2]. Thus, the analysis of the superhorizon evolution of the curvature perturbations is 
extremely important. 

Conventionally, the standard perturbative approach has been adopted to investigate 
the curvature perturbations [3]. Linear analysis gives simple equations and solutions, which 
are used to evaluate the power spectrum of the curvature perturbations. However, recent 
observations are precise enough to probe the nonlinear effects such as the bispectrum and 
the trispectrum of the curvature perturbations, which are crucially important to identify 
the origin of the primordial density fluctuations. Though all of the single-field inflation 
models with the canonical kinetic term predict negligible nonlinear effects on the curva- 
ture perturbations [I] (see also [5] for a single-field high friction model), inflation models 
with noncanonical kinetic terms [6] and light field models [7] such as curvaton [SJ E] and 
modulated reheating [TQl [11] can predict large non-Gaussianity of the curvature perturba- 
tions. Therefore, the standard perturbative approach must be extended to deal with the 
nonlinear effects. But the equations and their solutions often become too complicated to 
be understood intuitively. Then, instead of the conventional perturbative approach, two 
different formalisms are proposed to deal with the fully nonlinear curvature perturbations. 
One is the <5N formalism [12] and the other is the covariant formalism [131 114 j^i 

The SN formalism provides a powerful and simple method to investigate the super- 
horizon evolution of the fully nonlinear curvature perturbations. According to the <5N 
formalism, the superhorizon curvature perturbation on a uniform energy density hyper- 
surface at a late time is equal to the perturbation in the time integral of the local expansion 
from an initial flat hypersurface to a final uniform energy density hypersurface. At leading 
order in gradient expansion, which corresponds to the separate Universe approach [16], 
each superhorizon sized region of the Universe evolves in time independently from other 
regions. Thus, we have only to evaluate the expansion of the unperturbed Friedmann 
Universe. In fact, it is quite easy, by the use of the <5N formalism, to show that the fully 
nonlinear curvature perturbation on a uniform energy density hypersurface is conserved 
on superhorizon scales if the pressure is only a function of the energy density |17j . 

On the other hand, the covariant formalism defines covariant quantities corresponding 
to physical quantities and derives their evolution equations. Since these quantities are 

* 1 The space gradient of the fully nonlinear curvature perturbations was first introduced in |15j which 
is based on the leading order of the coordinate-based gradient expansion. It was later extended to the 
curvature covector in a covariant manner by [131 114) . 
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tensor, they are independent of a particular selection of a coordinate system and easy to 
understand from a geometrical perspective. In particular, the fully nonlinear evolution 
equations of the curvature covectors corresponding to the curvature perturbations are 
easily derived. Notice that, different from the SN formalism, they are exact and valid at 
all scales. Its covector is shown to be conserved for adiabatic perturbations [T3] . 

Then, one may wonder what is the relation between the above two formalisms. In 
particular, what is the relation between the curvature perturbation and the curvature cov- 



ector? At the linear, second and third orders [131 (HI HH] cfl the curvature covector is shown 
to be related to the gradient of the curvature perturbations of the standard perturbation 
theory [20]. For the case of a barotropic fluid, the relation at the fully nonlinear order has 
been discussed in [21]. One of the present author (Y.W.) has evaluated the bispectrum of 
the curvature perturbations in two-field inflation models by two different formalisms and 
explicitly shown that both of them coincide [22] . Similar analyses in ekpyrotic models for 
the bispectrum [23] and the trispectrum [TJ] have been done, but they have less accurate 
quantitative agreement than [22]. Nevertheless, to some extent, the SN formalism and the 
covariant formalism is shown to be equivalent for superhorizon curvature perturbations 
up to third orders. In this paper, we explicitly show the fully nonlinear equivalence of the 
SN and the covariant formalisms for the superhorizon curvature perturbations. 

The organization of this paper is as follows. In the next two sections, we briefly review 
the SN and the covariant formalisms, and derive the evolution equations for the curvature 
perturbation and the curvature covector from the continuity equation. In Sec. HI we show 
the fully nonlinear equivalence of the two formalisms for the superhorizon curvature per- 
turbations. In Sec. |5j isocurvature perturbations in the covariant formalism are discussed. 
We clarify the relation between the fully nonlinear evolution of the curvature covector 
and that of the curvature perturbation for multiple interacting fluids. The final section is 
devoted to the summary. 



2 5N formalism 

In this section, we briefly review the fully nonlinear version of the SN formalism, picking 
up the relevant points derived in |17j . 

Let us start from the Arnowitt-Deser-Misner (ADM) decomposition of the metric, 

ds 2 = -AT 2 dt 2 + ^{dx 1 + (3 l dt)(dx 3 + j dt). (1) 

We further decompose 7^ as 

ll3 = a 2 (t)e^(e%, (2) 

where a(t) is the fiducial scale factor, hij is a traceless tensor, and ip is the perturbation 
of a unit spatial volume. 

# 2 In Ref. [18], they define a scalar quantity similar to the curvature perturbation £ and discuss its 
evolution at all orders. 
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The gradient expansion assumes that all the quantities of interest are smooth over very 
large scales [23]. Practically, this expansion is achieved by a prescription that a spatial 
derivative acting on any variable is accompanied by a small quantity e and we expand all 
the equations in terms of e. It leads to the so-called 5N formalism by keeping only terms 
at zeroth and first orders in e and dropping higher order terms. As in the literature, we 
assume that 

F = 0(e), ^ = C(e 2 ), (3) 

where the prime (') denotes differentiation with respect to the cosmic time t. 

Assuming that the anisotropic stress is 0(e 2 ), which is valid in many cases, we can 
write the large-scale energy-momentum tensor of matter filling up the Universe as 

T vv = (p + P)u^u u + Pg^ + 0(e 2 ). (4) 

The four- vector is normalized to be = — 1 and the integration curve along gives 
the fluid worldline. Expansion of u 11 along the worldline is given by 

Q = V^ = ^- r (- + A+0{e 2 ). (5) 



Correspondingly, the local expansion rate H measured by the local observer using his 
proper time is 

H = lG = ±-(-+iP')+0(6 2 ). (6) 



3 N 



a 



The e-folding number of the expansion rate along the worldline is defined as 



N(t 2 ,t 1 ;x t ) = - / dt ATQ = log + ipfax*) - i>(h,x l ) + 0(e 2 ). (7) 
Jti a {ti) 

Therefore, if the t = const, hypersurface at initial time t\ is taken to be a flat one, 
i. e. ip(ti,x l ) = 0, this formula states that perturbation of the e-folding number from t\ to 
t 2 coincides with the curvature perturbation at t = t 2 (up to first order in e). Equation 
(0) is the essence of the SN formalism. 

Time evolution of ip can be derived from the continuity equation m^V^T 1 ^ = 0; 

^1 + j/ (t J) = P^Ml + (e 2 ) (8) 

a(t) + W> X) 3[p(t,^) + P(t,^)] +U[€ h ^ 

For a fluid having a barotropic equation of state P = P(p), we can integrate this equation 
to obtain a conserved quantity £ given by 

1 do 
3 J P (t) P + P{p) 

In particular, if we take a uniform energy density hypersurface E p , the second term in 
the right hand side vanishes. Therefore, the curvature perturbation on that hypersurface 
remains constant in time. 



# 3 Mathematically speaking, the second assumption is not needed to derive the evolution equation of ifj 
in Eq. ([5]) but needed to interpret tp as the spatial curvature. 
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3 Covariant formalism 



In this section, we briefly review the covariant formalism developed in Refs. [H2 S3]- 

While the SN formalism takes the coordinate-based approach from the outset, the co- 
variant formalism first defines fully nonlinear and covariant quantities and then derives the 
basic equations. After this, we can choose a coordinate system and expand the quantities 
and equations up to the desired order. 

In the covariant formalism, the curvature perturbation is represented by a covector ( a 
defined by 

C = d a N - -d aP = d a N - ^d a p. (10) 
P 3p 

We use a, b, ■ ■ ■ to stress that we are working in the covariant approach. Here N is the 
local e-folding number along the fluid worldline introduced by Eq. (JJJ) up to an integration 
constant that vanishes on some hypersurface. In the covariant form, it is given by 

N =- fdrQ, (11) 



3 . 

where r is a proper time of the fluid. The continuity equation u a VftT b a = becomes 

p + 0(p + P) = O, (12) 

where we have assumed that the Universe is filled up with a perfect fluid. Using this 
equation, we can derive the evolution equation for ( a along the fluid world-line; 

^-wrp){ d ° p -j d *^ (13) 

Here the dot Q denotes the Lie derivative along the fluid worldline. For example, 

P = U a V a P, ( a = U%( a + ( b d a U b . (14) 

It is worth stressing that the evolution equation (fT3j) is exact and valid on all length scales. 
For fluids having a barotropic equation of state P = P(p), the right-hand side of Eq. (fT3|) 
vanishes and ( a is therefore conserved along the worldline. 



4 Equivalence between the two approaches 

Now, let us again choose a coordinate system given by Eq. flTJ with the condition (j3J) and 
expand ( a on this coordinate. For the time component, we find 

^ N ' + m0k^r o{c2) ' (15) 



l 



where we have used the continuity equation (fl2"|) to show the last equality. For the spatial 
component, we have 

& = diN+ i\ 77 J^pf. K] 9ip(t,x j ). (16) 
3[p(t,x J ) + P(t,x t )\ 



Although the right hand side is not generally a total derivative^ 4 !, it becomes a total 
derivative on a uniform energy density hypersurface £ p . Therefore, we find that Q is the 
spatial derivative of the perturbation of the e- folding number evaluated on E p . Using 
Eq. (GO), Q can be written as 

C<(*2, x') = <W(t 2 , x j ) - tP(t u xi)) + 0(e 3 ), (17) 

where ^(t 2 ,x- 7 ) must be understood as the curvature perturbation on S p . In particular, if 
we take the flat hypersurface at t — t\, the above equation reduces to 

(,fe^) = ^(i 2 ,^) + 0(£ 3 ). (18) 

At this stage, it is clear that Q on very large scales is merely a spatial derivative of the 
curvature perturbation on £ p . 

We can explicitly show that the evolution equation (TT3T) evaluated on S p indeed leads 
to Eq. ([8]). The left hand side of Eq. (TL3]) is given by 

C i = ^d^ + 0(e 3 ). (19) 
On the other hand, the right hand side is given by 

' d t p- -d lP ) = -—EMl-d.Pit,^). (20) 



3(p + P) \ p ' ) p(t)+P(t,x 

Then using a formula for M valid on E p , 



M = P^ +P ^ (21) 

p(t) + p(t,x*y {2l} 



# 4 One exception is for the fluid having a barotropic equation of state P = P(p). In [21] it has been 
shown that 

logp 



& = d t 



N ■ 



3(1 + «>). 

where w = P{p)/p is constant in space. Integrating over space and choosing the integration constant 
properly, one can immediately get 

C = 5N+\ f ^ P 5N = N-N(t), 
3 Jp{t) (1 + w)p 

which is equivalent to Eq. ©. In the following, wc do not demand the barotropic condition. 



and a fact that S p and the uniform Hubble hypersurface coincide to first order in e [TT] . 
we find that Eq. ( !T3|) reduces to 



di4)' = di H{t) 



p(t) + P(t) 



(22) 



p(t) + P(t,xi) 



whose integration over x l , combined with the continuity equation p'(t) /M = —3H(t) [p(t) + 
P(t,x 1 )], coincides with Eq. (jSJ provided the integration constant is chosen properly, i.e., 
a! I a. The curvature perturbation associated with its covector obeys the same evolution 
equations as for the one given in the <5N formalism, which therefore explicitly establishes 
the equivalence at the fully nonlinear level between the 5N formalism and the covariant 
formalism. 

5 Isocurvature perturbation in the covariant formal- 
ism 

Finally we will show that the evolution equation (fl3|) . if the fluid consists of multiple 
interacting fluids, can be written in such a way that the right-hand side becomes the sum 
over all the possible combination of the isocurvature perturbations between two different 
components, which resembles the standard one obtained under the linear approximation. 
By setting a coordinate system, we clarify the relation between the isocurvature covector 
and the isocurvature perturbation at the fully nonlinear level. 

To rewrite the right hand side of Eq. ( lT3"j) in the desired way, we first introduce the 
curvature covector for the fluid component A by 




(23) 




(24) 



We can then rewrite the right-hand side of Eq. ffTB"]) as 



3(P + P) 







d a P - —d a p 



p 



) 



3(P + P) 



e 




p A P B S A a B X2h) 



where we have assumed that there is no dissipation in total and 




(26) 



6 



is the sum of the intrinsic nonadiabatic perturbations for each fluid. Therefore, Eq. ( [13]) 
becomes 

This result is a fully nonlinear generalization of [26] and obtained in Ref. [25]. One can 
clearly see that the evolution of the curvature covector is generated by the existence of 
isocurvature covectors and the relative difference of the sound speeds of fluids. 

It is important to notice that we can apply this equation even when energy is exchanged 
among the fluids. To see it in a more manifest way, let us introduce the effect of the energy 
transfer by 

p A + e(p A + P A ) = Q A , (28) 

where Q A represents a rate of energy density that the fluid A absorbs. Since the energy 
is conserved in total, we have ^ A Q A = 0. Using this equation, we find that Eq. (|27|) can 
be written as 

• _Q^fP A _P B \p A + P A p B + P B { d a p A d a p B 



6 V P A P B P + p P + p \P A + P A p B + P B 



AB 

Q2 



3p2 



- (pT^ + (Q)d a p , (29) 



where (Q) = Y1a c aQ A - As * s clear, the last term represents a contribution from the 
energy transfer but the energy transfer also affects the other terms implicitly through the 
evolution equation (|2"8"|) . 

If each fluid obeys a barotropic equation of state [Fa = 0] , by applying the coordinate 
system given by Eq. ([1]) to Eq. ( 129|) . the evolution equation of the curvature perturbation 
corresponding to Eq. (129]) can be derived. Choosing again the uniform energy density 
hypersurface, the last term of Eq. ( 129]) vanishes and we end up with the following equation: 

6 ^ p+P p+P 

where S AB = 3(( A — £ s ) is the fully nonlinear isocurvature perturbation introduced in 
and C A is defined by 

■An jt\ - ,.u J, . 1 f A( '"' ) ip' 



ew^sn-j^ (3D 

Thus, we have shown that the evolution of the curvature covector in the covariant formal- 
ism is equivalent to that of the standard curvature perturbation for multiple interacting 
fluids as well. Notice that P A /p A coincides with the sound velocity squared c\ for a 
barotropic fluid. 
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6 Summary 



In this paper, we have shown the fully nonlinear equivalence of curvature perturbations 
on superhorizon scales from the two formalisms: the SN and the covariant formalisms. In 
particular, by setting a coordinate system and integrating over space, we have identified the 
evolution equation of the curvature covector with that of the curvature perturbation in the 
5N formalism. The key assumption here is that the matter energy-momentum tensor takes 
the perfect fluid form on superhorizon scales, i.e., the continuity equation ( TT2|) completely 
describes the nonlinear evolution of the energy density of the system. 

We have also clarified the relation between the isocurvature covector and the nonlinear 
isocurvature perturbation for multiple interacting fluids. This identification enables us 
to bridge the results of standard perturbation theory to the covariant formalism. Our 
treatment here is somewhat less general than (25] but it is enough to obtain the important 
relation to the standard curvature perturbation. It would be interesting to extend our 
treatment to systems with multiple scalar, spinor and vector fields, where the next order 
in gradient expansion may become important. 
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